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We report on a first, comprehensive partially quenched study of the r\ — rj problem, based on SESAM config- 
uration on a 16'^ X 32 lattice at /3 = 5.6 QCD with two (mass degenerate) active sea quark flavours, N j — 2. By 
means of the spectral approximation of the two-loop (hairpin) diagrams, we find clear plateau formation in the 
effective masses which enables us both to determine the -q-ri' mass matrix and the a-parameter in the effective 
chiral Lagrangian for the fiavour singlet sector, a — 0.028 ± 0.013. 



1. INTRODUCTION 

Flavour singlet pseudoscalar states bear wit- 
ness to the topological structure of the QCD vac- 
uum and hence constitute objects of great physi- 
cal interest. In lattice gauge theory, considerable 
progress has been achieved recently 'j^ in directly 
accessing such 77' like states from their propaga- 
tors in Nf = 2 QCD with Wilson fermion ac- 
tion, based on applying spectral approximations 
to two-loop correlators (hairpin diagrams). 

In this note, we extend this work to the more 
realistic setting of partially quenched QCD with 
two active sea quark flavours, where valence and 
sea quarks are allowed to carry different masses. 
In our comprehensive study we analyzed some 800 
independent SESAM configurations on 16^ x 32 
lattices at /3 = 5.6 P], with four different sea 
quark masses, as indicated in Fig^^. 

2. The partially quenched scenario 

Partial quenching in the flavour singlet sector 
of QCD is described in the framework of chiral 
perturbation theory by a nonsymmetric piece of 
the effective chiral Lagrangian PJ 
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Figure 1. Our points of operation in the hopping 
parameter plane of valence/sea quarks. 



This expression introduces two effective 'low 
energy constants': the mass gap ^ and an inter- 
action parameter accounting for momentum de- 
pendence, a. The Lagrangian translates into two- 
point correlation functions (propagators) for the 
neutral pseudoscalar quark bilinears 
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where the index i runs over the quark and pseu- 
doquarks d.o.f., the latter being geared to act as 
'determinant eaters' on the valence quarks 

In momentum space the disconnected part of 
the ps-correlators has a compact form |^ with 
the infamous double pole structure 



'^For further details of tfie analysis, see Refs. I1I5I . 
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Figure 2. First evidence of plateau formation in 
a partially quenched setting, illustrated on the 
effective mass gap, morm(i), in the disconnected 
two loop correlator, Z?„„(i). The data refer to 
the lightest sea quark mass of SESAM, i.e. Kgea = 
.1575, for various isosinglet valence quark masses. 



Figure 3. Extrapolations of the quadratic mass 
matrix elements of light (n) and strange (s) isos- 
inglet valence quarks to QCD with chiral sea 
quarks (Eq. QIl, at a = 0. 



The prefactor, Af = 1/(1 + Nfa) and the singlet 
mass, M^, = (Mj + NfiJ?)/i\ + Nja) carry ex- 
plicit a- dependencies. Rememeber that the octet 
masses {Mii,Mjj, and Md) from valence and sea 
quarks, respectively, are readily computed from 
connected diagrams on the lattice. Our present 
concern is the determination of the mass gap, /i, 
and a from our lattice data. 



3. Results 

As it has not been possible to determine the 
effective constant a in the past we explore the 
potential of our data for achieving progress on 
this issue by applying our spectral techniques ^ 
in two steps. 

3.1. a = analysis 

In a first attempt, we choose a = and study 
effective mass gaps as determined from the hair- 
pin diagrams, Dij{t), by fitting, w.r.t. raQij{t) := 
fi at each value of t, the zero-momentum Fourier 
transforms of ansatz Eq. |21 to our data. In bro- 
ken SU(3) with two active quark flavours, we in- 
troduce two isosinglet pseudoscalars, designated 
by the indices n (for 'light', nonstrange) and s 



(for strange). In this manner, we can study three 
types of hairpin diagrams, Dnn , Dns , Dss, with 
the strange quark mass determined on the lattice 
from the kaon mass (for k locations, see Fig.l). 
As a result, we find satisfying plateau formations 
in ■moij{t) as illustrated in Fig. [21 In fact, the 
numerical quality of our signals comfortably al- 
lows for the consecutive extrapolations (i) to light 
valence and then (ii) to light sea quark masses; 
the chiral sea quark extrapolation is exhibited in 
Fig. El 

From the consecutive chiral extrapolations, we 
obtain for the quadratic mass matrix in the 
quark-flavour basis 



M := 



O.sn 



nn nn ' 

mass scale 2 we find 



'Onn ■ 
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With the p-meson 



/ (306 ±91)2 ^/2(201±34)M 
V2(201 ± 34)2 (ggo ± 15)2 ) iVieV 

Diagonalization of M. renders 

= 292 ± 31 MeV, M^, = 686 ± 31 MeV. (5) 

These numbers from N f ~ 2 look promising: the 
individual mass values - while being low w.r.t. 
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Figure 4. Minimizing the valence quark depen- 
dence of the mass gap by variation of a. 



Figure 5. Univerval plateau formation, at the 
lightest sea quak mass, adjusting to a = 0.0328. 



the real three flavour world - yield a mass split- 
ting that compares very well to phenomenology. 

3.2. First determination of a 

We remark that the outcome of the analysis at 
this stage is not in accord with the tree approx- 
imation to the Lagrangian, Eq. |2 as the latter 
predicts the mass plateaus to be independent of 
the valence quark mass, m^,, contrary to the find- 
ings in Fig. 121 Therefore, we repeated the above 
mass plateau study on a set of nonvanishing a- 
parameters, in an attempt to verify the compli- 
ance of our data with such independence: inspect- 
ing the plot of ruonn vs. m„ (Fig. 0J it becomes 
obvious (i) that TOorm(mt,) shows a simple, linear 
behaviour and (ii) that a can indeed be adjusted 
to produce a zero slope of monni'mv) ^! 

In this way it is straightforward to find an op- 
timal value, aopt, that eliminates (for a given sea 
quark mass) the m.u-dependence from monn, re- 
sulting in a universal plateau level. This collapse 
of data into universal plateaus (of heights jl) is ex- 
emplified in Fig.|S| The comparision with the sit- 
uation encountered with a — (Fig. ^ provides 
clear evidence for our sensitivity in determining 
e^opt • 

After chiral sea quark extrapolation we arrive 



at the values 

a == 0.028 ±0.013 



^ = 203 ± 34MeV 



(6) 



^The situation is less favourable for the gap values ex- 
tracted from Dn«l51. 



Note that the numerical value for the mass gap 
at the chiral point is robust w.r.t. to the two dif- 
ferent approaches presented here. One might say 
that the allowance for the a parameter anticipates 
the chiral valence quark limits. 

In conclusion: the application of spectral meth- 
ods to the hairpin diagrams allows rather detailed 
partially quenched studies in the flavour singlet 
sector, with standard Wilson fcrmions. This pro- 
vides a strong motivation to extend this work to 
the case of overlap fermions. 
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